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Abstract 

We studied the high energy elastic scattering of quark anti-quark with an 
exchange of a mesonic state in the t channel with — i/A 2 S> 1. Both the 
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■ normalization factor and the Regge trajectory can be calculated in PQCD in 

cases of fixed (non-running) and running coupling constant. The dependence 
of the Regge trajectory on the coupling constant is highly non- linear and the 
trajectory is of order of 0.2 in the interesting physical range. 

(Submitted to Physical Review D) 



1 Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



1 Introduction 

With the advances of LHC and SSC, it is possible to study experimentally the Regge 
behavior in the parton level where the momentum transfer squared —t ^> A 2 but is 
still smaller than the center of mass energy squared s of the partons. The Regge 
limit of the parton scattering amplitudes corresponds to the small x limit of parton 
distribution, while the presence of a "large" scale —t justifies the use of perturbative 
QCD. 

The sea-quark and gluon distribution for small x is related to the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) Pomeron The distributions of sea-quarks and gluons 
grow like x~ ap at small x. where ap is the trajectory of the BFKL Pomeron. The 
behavior of the valence-quark distributions is controlled by the mesonic Reggeons Q. 
It grows as x~ WR with ujr the trajectories of the mesonic Reggeons. However, great 
care is necessary to separate the perturbative behavior from the non-perturbative soft 
physics. 

Several hard partonic processes to measure the behavior of the BFKL Pomeron 
[2-7] have been discussed in the literature. In Ref f7j, A.H. Mueller and the author 
proposed a process of high energy, fixed t parton-parton scattering through the ex- 
change of a BFKL Pomeron. It is natural to extend this idea to the process whereby a 
mesonic Reggeon is exchanged. It is the objective of this paper to set up the necessary 
machinery to investigate the process of the mesonic Reggeon exchange and to study 
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the possibilty of measuring such behavior. The scattering amplitudes are calculated, 
while the normalization factor and the trajectory can be obtained explicitly. 

In this paper, we study the quark anti-quark scattering amplitudes [Fig.l] with 
flavor exchange in the t channel. In the kinematic region, 



where the auxiliary parameter /i is the infrared cutoff of the transverse component of 
the momenta (with respect to the initial momenta pi and P2) of virtual particles in the 
Feynman integrals where the terms ~ a(/i 2 )[(a(/z 2 )/7r) ln 2 (s//x 2 )] n are summed, and 
a(fi 2 ) is the strong coupling constant at the scale /1 2 . This is the double logarithmic 
(DL) approximation. The method of separating the softest virtual particle || allows 
one to calculate the partial wave amplitudes in the double logarithmic approximation. 
In Ref ||, R. Kirschner and L.N. Lipatov give equations for the partial wave of the 
amplitudes for both the color singlet and the octet exchange in the t channel. The 
octet exchange is suppressed, because it has a strong tendency to radiate gluons in 
the rapidity interval defined by the colliding quark anti-quark. This phenomenon is 
reflected by the negative intercept of the Regge trajectories of the octet exchange. In 
what follows, we restrict ourselves to the consideration of the color singlet exchange. 
Based on the partial wave results in Ref ||, we study in detail the scattering am- 
plitudes for cases of both fixed (non-running) and running coupling constants, using 
analytic methods as far as possible. 



s ~\ u |> -t ~ /j 2 > A 2 , 
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The outline of the paper is as follows: In Sec. 2, we review the results of the 
partial wave amplitudes for the color singlet exchange. Scattering amplitudes for the 
fixed coupling case are presented in Sec. 3 for both positive and negative signature 
cases. In Sec. 4, both approximate and numerical methods are employed to study the 
effect of the inclusion of a running constant in the positive signature. Finally, Sec. 5 
summarizes our conclusion. 



2 Review of the partial wave amplitudes 

We consider the amplitudes of annihilation (qq — > QQ) with the exchange of a meson- 
like state in the t channel. With respect to the color group SU(N), the amplitudes 
can be decomposed into a singlet [M (s)] and an octet [M„(s)]. In double-logarithmic 
approximation, the spinor structure of the Born term is preserved in higher order, so 
we can write the amplitudes as b M(s, t) where b = 7^ <8> 7^ / s is the Born amplitude 

but without the coupling constant g 2 . M Q and M v are 

N 2 — 1 2 1 2 

M \ B orn= ^ 9 , M v \ B orn= ~~^jy9 , ( 2 ) 

after using the color projectors P (singlet) and P v (octet) 

Polb = -jy^aa'Sw , Klb = ^ab^a'b' ~ Saa'tibV (3) 

and the decomposition 

nfa'b' _ p a'b' yr , p a'b' yr I a\ 

m ab — r Q ab 1VI + r v ab lvl v K^J 
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where a, b and a', b' label the color states of the initial and final quarks. 

The asymptotics of the scattering amplitudes at large s and fixed t are determined 
by the singularities of the partial wave fj(t) in the crossed channel. In order to express 
the amplitude in terms of partial wave, we need to divide the amplitude M into 
parts that are symmetrical and anti-symmetrical with respect to the transformation 
s <-> u ~ —s : 



In double-logarithmic approximation, the Sommerfeld- Watson transformation reduces 
to the Mellin transformation, and the even (odd) part of the amplitude is related to 
the positive (negative signature) partial wave 



with f p {uj) includes the factor (sinirui) 1 usually written explicity in the Sommerfeld- 
Watson integral. The signature factor is given by 



Using the method of isolating the softest virtual particle with the lowest transverse 
momentum k\ in the Feynman diagrams, R. Kirschner and L.N. Lipatov were able 
to give equations for the partial wave amplitudes. The positive signature amplitudes 




(5) 
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are 



while the negative signature amplitude is more complicated, 
with 

TV 2 — 1 1 
a ° = 2N ' a " = ~2iV' 6 " = ^ ^ 

and boundary conditions 

tfUx, = ^ (< = 0, W ). (11) 

Here we list the partial wave amplitudes that are relevant to the color singlet 
exchange. The first terms in the above equations are the contribution from the Born 
terms which have pole singularity at u = 0. 

The equations of color singlet exchange are purely algebraic, so that the solutions 
can be written in explicit form: 

tf(u) = 4^ [ 1 - y/1 - (£) 2 ] 

/ » = 4^ [ 1 - - (£) 2 (1 - ^/+) ] (12) 
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which shows that fo~{u)) has a square- root branch point at (Fig. 2) 

while fo(oo) has singularity to the right of c<j + so that it becomes dominant when 
energy is large. The equation for f^{uj) is of Riccatti type. It can be solved easily 
and leads to the following result: 

/+(.) = ;v£ln[e*P(i(^>,(-)] 

* = * = ^ = -2^- (14) 



where T> u (z) is the parabolic-cylinder function JIT . 

In deriving fi(oo), the coupling constant was taken to be fixed. But as the trans- 
verse momentum k\ covers the large range from /i 2 to s 0, contribution from the 
running coupling is important and cannot be neglected. As demonstrated later, the 
inclusion of the running coupling constant changes the singularity of the partial wave 
amplitude from square root singularity to pole singularity. 

Taking into account the effect of running coupling, the positive-signature singlet 
channel partial wave amplitude F{oj,h 2 ) satisfies 

^^Fiu^^-uF^^ + ^Fiu^^f + aog 2 ^ 2 ) = (15) 
with the boundary conditions 



F(uj,h z ) -t . (16) 



The solution to the above equation is given by 



,2, _ ,2 ^ i_r,T,/ fl n . W 



F(^) = 8.V^ln[*(-3^,0;^)] (17) 

where b = (yiV— 2/3iV/)/167r 2 , the coefficient of the first term of (3 function, \I/(a, c; z) 
is a confluent hypergeometric function [[[(]] defined by 

q(a,c,z) = — - / e-**t a - 1 (l + t) c - a - 1 dfc. (18) 
r (a) jo 



3 Fixed coupling 

For the positive signature channel, 



ra+ioo A, , — o 
J a—ioQ Z7T 6 6 



where we have set /i 2 = — t. Substituting fo(uj) from eq. (|I2"D, we have 

ra+ioo A, , — c , ■ 

M + (s) = / -^(-^4^ [1 - - (^)% (20) 

where the contour of integration is shown in Fig. 2. We take the branch cut from — uo + 
to so the contour of integration can be deformed to encircle the branch cut. 

The positive signature amplitude becomes 

M + ( S ) = ^ (^r 4^ Dis [1 - v/l-(^) 2 ] • (21) 



With Dis [1 — Jl — (^) 2 ] = 2i ^ + ^ I 1 and modified Bessel function I v defined 
by 

/ (V) = (*/ 2 )" f 1 (i - t 2 Y- l / 2 e±- zt dt (22) 

A ) r(z/ + i/2)r(i/2)7-i l ] 6 d [ ] 

Mq(s) is evaluated as 

M+(s) = ^(w+y) (23) 

where ?/ = ln(s/ — t), the rapidity interval between the quark anti-quark pair. If y is 
large, i.e. in the asymptotic region, 

h{u) + y) ~ 



\/27ru; + 7/ 
and 

which shows the Regge limit behavior. 

As mentioned before, the negative signature channel dominates asymptotically, so 
an explicit solution is desired. However, because /+ depends on the parabolic-cylinder 
function, it is not possible to have a convenient solution. In view of that, we take 
lim N — > oo. The limit is not just an academic exercise but has its own significance, 
because the relevant parameter in this problem is p v = —1/2N 2 = —1/18 which can 



be taken to be zero without introducing much error. With [11 



dz 



fy(uj) can be rewritten as 



fAu) ~ ~^-^m~- (25) 



In the limit, p v = -1/2N 2 -> while 



^ 2 = a - (iV2 - 1) ^^ = 2^. (26) 
7riV 7r 



With the help of 0, 



V {z) = e~ z2/ \ and 
T>-i(z) = y^e^ 4 Erfc(^), (27) 

where Erfc(z) is the complementary error function. can be simplified as 

/+(-) = e<4> a Erfc(£). (28) 

Define u = — , and drop the first term of the negative partial wave amplitude 



in eq.( 0) as it has no singularity and hence contributes nothing to the amplitudes 



after performing the Mellin transformation. f (u) then reduces to 



----- LJ y,<r>- ii - -^-j < "'V.vUiu) . (29) 



u 



Expression (29) is corrected up to the order 1/N 2 and clearly the correction term 
comes from f+{u). The zeros of the expression u 3 — u — e u2 Erfc(u) evaluated 
numerically are found to be u + = 1.0388 and m_ = —0.4668. As u + > 1, the square 
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root singularity of the negative channel lies to the right of that of the positive signature 
channel. But u + — 1 = .0388 is not large. This indicates that both channels have 
similar behavior phenomenologically. The expression u 3 — u — 2 ^-e" 2 Erfc(M) can be 
approximated by (u — u + ){u — u ) 2 within the region u_ < u < u + which is chosen 
to be the branch cut. With this approximation, 



(2^+ r~ 

u v 



(30) 



Performing the Mellin tranformation by deforming the contour of integration to 
enclose the branch cut from u = to u = u + (Fig. 3) leads to 

Mn(s) = 7t 2 (cu+) 3 /" + du (—Y +u ( u - u-) BisJu(u - u+) 
Jo t v 

= 2i7r 2 (cu + ) 3 du (— T u (u-uJ) | u(u-u+) | 1/2 . (31) 
Jo t 

Change variable v! to u — u+/2, and the integral can be expressed in terms of the 
modified Bessel function I u , 

Mo(s) = I'—e^l 2 [(^±-u--^)I 1 (u + u + y/2) 
V 2 io+y 

+ ^±I ( u+ cu + y/2)}. (32) 
In the asymptotic region where y is large, the amplitude reduces to 

M- y (s) = l ^{u + - «_) u + e^- 1 ^ M+ y {s) (33) 



with M+ y (s) given by eq. (|4j). In the range of SSC energies where y ~ 8 and 
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uo+ ~ 0.5 



^{u+ - u-) uj + e {u+ ~ 1)uJ+y ~ 1, (34) 

which suggests that the positive and negative signature amplitudes are of the same 
order and cannot be distinguished in the interesting energies range. Actually, y needs 
to be about 50 before the negative signature channel is appreciatively different from 
the positive signature channel. The difference between y/2u v and uo + in eq. (|26 ) 



enhances the negative signature channel a little bit. Effectively, it changes u + — 1 = 
0.0388 to 0.1, but its effect is still small at SSC energies. 

The two amplitudes have a phase difference of ~ ir/2, because M+ y (s) is purely 
real while M~ (s) is purely imaginary, and they have nearly the same magnitude. 



Both give a non-linear quark anti-quark Regge trajectory u(t) ~ Ja s (—t). 



4 Running Coupling 

As we have already mentioned, there is reason to believe that the inclusion of running 
coupling effects will greatly change the behavior of the amplitudes. In this section, we 
would like to restrict ourselves to the study of the positive signature singlet channel. 
It is plausible to suggest that the inclusion of the running coupling effect would not 
change the similarity between the positive and negative signature amplitudes greatly 
for the following reasons: the only difference between the positive and negative chan- 
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nels stems from the contribution of the double logarithmic soft gluon; the soft gluon 
contributes to the negative signature singlet channel but not the positive signature 
singlet channel ||. However, as seen in the previous section, the soft gluon does 
not change the amplitude appreciatively. The inclusion of the running coupling will 
decrease the importance of the soft gluon contribution as the coupling between quark 
and gluon is smaller than in the fixed coupling case. Therefore, for practical purposes, 
we can take positive and negative singlet amplitudes to be equal in magnitude but 
with phase difference tt/2. 

In the following subsections, we will study the effect of the running coupling in 
the positive signature channel using two different methods. In the first method, we 
take lim b — > 0, where b is the coefficient of the first term in the (3 function. In 
this approximation, we will recover the previous fo~(uj) result with a correction term 
which is liner in b. However, as we will show later, the correction term enhances 
the amplitude by (lny) 3//2 relative to the fd(uj) term, so that the correction term 
dominates asymptotically and the approximation breaks down. The breakdown of 
b expansion in the asymptotic region closely relates to the fact that the transverse 
momentum in the loop integral extends to the order of s. It illustrates that in the 
asymptotic region, the running coupling effect is not a small perturbation to the fixed 
coupling, but rather, it changes the behavior of the amplitudes dramaically. 

In the second method, we take the approximation A = ao/(87r 2 buj) — > oo. It gives 
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a fairly accurate result with ~ 10% error compared to the numerical calculation, even 
though the value of A lies between 1 and 2 when the running coupling constant is in 
the interesting physical range. 
Let us define 

Both A and v are positive if u > 0. In the color singlet channel, the trajectories lie 
to the right of the origin of the complex uj plane, because the higher order diagrams 
enhance the amplitudes. This is in contrast to the case of the color octet exchange 
which suppresses the amplitude. It is safe to assume that both A and v are positive. 
As ^'s dependence on /i 2 is through the running coupling constant # 2 (/^ 2 ), it is better 
to write 

d 2 dg 2 d 



<9/i 2 d/j, 2 dv 
where /-t 2 f^2 = — bg 4 . In terms of A and z/, the partial wave amplitude F(uj,j2 2 ) is 

F(u, jj 2 ) = 87T 2 w^ln[^(-A,0;z/)]. (37) 



4.1 b expansion 

In the limit b — > 0, A, v — > oo ; ^ can be approximated by [|T0 



tf(-A,0;i/) = AV 1/4 (z/-4A)- 1/4 

exp(-A + \ - - AXfl 2 + A In [1^^]) (3? 
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which leads to 

|-M.(-a,M] = ^-d-v' 1/2 -I^4A)l- < 39 > 

The second term in the above expression turns out to be independent of b and with 
the first term, it recovers f^(u). The third term, which is the correction term, is 
found to be 

2 A 2nba s Lj +2 

It is linear in b and has a pole singularity at u = u + in contrast to the square root 
singuarity of the first term. Putting eqs. fl3~T|), (p9[) , and (^) together, we get 

FH = ,41) 

The negative sign of the correction terms indicates that the running coupling reduces 
the overall amplitude. This behavior is also found in the Pomeron exchange [|J. The 
contribution of the correction term to the positive signature amplitude is 

AM+(s) = -^-ba s u + e» + v (42) 

which leads, in the asymptotic region, to 



M+ n = M+ + AM+ 

[1 - {2^\ba s u + y^}Mt sy (43) 
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where the correction term has y 3 ^ 2 dependence relative to the first order term. The 
approximation is valid only when 

(2n) 1/2 irba s LU + y 3/2 < 1. 

With uj + ~ 0.4, a s ~ 0.20 and b = 0.053 (four flavors), the above inequality implies 
y <C 9.7 which is the boundary of the asymptotic region. 

4.2 A — > oo approximation 

Before taking any approximation, let us study the nature of the singularity of the 
partial wave amplitude F(uj,fi 2 ), given by eq. fl3~T|). Both A and v are positive pro- 
vided that u, the Regge trajectory, is positive. Using the differential property of the 
confluent hypergeometric function [111, 



d 

— #(a,c, ;u) = -a*(a,c+l;i/) (44) 



F(u, fi 2 ) becomes 



2 A*(l-A,l;i/) 



F <^ 2 > = 8 ^ wn^r (45 > 



^(a, c; z/) is a many- valued function of u, and we usually take its principal branch in 
the plane cut along the negative real axis. Therefore, ^(a, c; v) is analytic for v > 0. 
For v > 0, the singularity of F(u,fi 2 ) must be the zeros of A,0;z/). \P(a, c; i/) 



cannot have positive zeros ifa>0orl+a — c>0 |]TD . This implies that 1 < A must 
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be true for A, 0; i>) to have positive zeros. From the definition of A, A ~ 1/u;, 



where 0J max corresponds to A = 1, and \1/ has one zero. Actually, according to Ref |TD 



the range 1 < A < 2, A, 0; u) has one zero. This can be seen by noting that |10 

q?(-n,c;u) = n\(-l) n L c ~ l (u) (46) 
where L^iv) is the generalized Laguerre polynomial [II]. Therefore 



tt(-l,0;i/) = (-1)^7 V) 



(47) 



as L~ n (v) = (—u) n /n\. The zero of \J r (— 1,0; z/) is at z/ = which in turn gives 
a s ^ oo as w = w ma3; = ao/(87r 2 6) = 0.32 is fixed by A = 1. For A = 2, which means 
w = u max /2 = 0.16, 

*(-2,0;i/) = -i/(2-i/) (48) 

where L^iy) = —v(2 — ^)/2 is used. Both u = and z/ = 2 are the zeros of 
^(—2,0;^), z/ = 2 corresponding to a s = u; max /(167r6) = 0.12. The v = solution 
is the subleading trajectory and so we do not consider here. Hence, from A = 1 to 
A = 2, a s ranges from 0.12 to oo which covers a large range of the interesting physical 
region. 

From the above discussion, we see that interesting physics already lies inside the 
range 1 < A < 2. A is of order 1, so it is not very large. Neverthless, let us take a 
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large A limit and compare the u(a s ) obtained with the exact results at a s = oo and 
a, = 0.12. 



When A — > oo, A, 0; u) can be approximated as JlO 



*(-A,0;i/) ~ 2 1 / 2 A A - 1 / 4 z/- 1 / 4 e^ 2 - A cos(A7r-2(Az/) 1 / 2 -^), (49) 

which in turn gives 

F(c,/i 2 ) = 8 7r 2 ^[-i ; + ^ + (^) 1 / 2 tan(A7r-2(Az/) 1 / 2 -i)]. (50) 

The first two terms can be dropped as they contribute nothing to the scattering 
amplitude after performing the Mellin transformation. From Eq. (OT), we see that 
F(u,f/?) has pole singularities at 

Avr - 2(A^ 2 - ~tt = ^vr, (51) 

where the leading trajectory corresponds to n = 1, 



r3 _|_ 2 / u) max \ 1/2] 
L 4 ' 7T V 47rfe« s / J 

~ r i + e^+i ■ ^) 

When a s = oo, u) Tun = 4u max /3 which is 33% larger than the exact value. For 
Oi — ^max/ (167T&) = 0.12, u run = 0.99uj max /2, which is only 1% off. The above analysis 
indicates that the ui run obtained from the large A approximation is accurate in the 
region where the running coupling is small. Equation (]52f) shows that the trajectory of 



the mesonic Reggeon is highly non-linear in a s , and one can see that u run (a s ) is always 
less than u + (a s ). The singularity of the partial wave amplitude in running coupling 
moves to the left of that of fixed coupling. For a s = 0.20, a typical value, u run = 
0.19 which is less than half of u + = 0.40. The smallness of the mesonic Reggeon 
trajectory imposes a serious problem of observing the mesonic exchange in the parton 
level experimentally, especially given that the mesonic exchange amplitude is already 
suppressed by a factor of s relative to the Pomeron exchange which constitutes a 
strong background. 

The amplitude can be obtained by the Mellin transforming Eq. (|50|), 



There is no justification for the correctness of the above expression in large A ap- 
proximation, although the trajectory does give fairly accurate results. Therefore, it is 
natural to consider the numerical evaluation of the amplitude and to compare it with 
the approximated result obtained in this section. To our surprise, the formula (^) is 
accurate to within 10%, as will be shown in the next section. 




(53) 
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4.3 Numerical Calculation 

Let us begin the numerical calculation by performing the Mellin transformation on 
F(u,fi 2 ) using the expression (flop . This leads to 



, dui ,—s.,, 9 A\I/(1 — A, 1: v) 



2th t *(-A,0;i/) 



\<k 2 lo (-^) w n(a a ) ^ (54) 



with 



n(a s ) = 



i _| g 2 b dty(\,0;v) \u=u R 

1 ' d(-A) 



(55) 

where co>#, the leading trajectory, is the root of \?(— A, 0; v). We are interested in the 
range 1 < A < 2. The roots un(a s ) in this range are shown in Fig. 4. The large A 
approximation suggests that we use the following parameterization to fit the data in 
Fig. 4, where a and j3 are constants: 

UR{a s ) = _ 1/2 (56) 

a + pa s 

Numerical fitting gives 

a = 1.14 x-, and (3 = 0.90 x 1(^)V2 (57) 

4 7T 47TO 

when we use the data from 1.5 < A < 2 corresponding to 0.121 < a s < 0.382. 
Compared with a = 3/4 and (3 = ^( i£ ^ L ) 1 ^ 2 in large A approximation, it once again 
confirms that it is a good approximation for the trajectory lor. 
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For the normalization n(a s ), the results for a s = 0.121 to 0.382 are shown in 
Fig. 5. They can be summarized by the following formula: 

n(a s ) = 0.8958tt6^(^^) 1/2 , (58) 

27rao 

where uji(a s ) is given by Eq. fl56|). Therefore, 

M r + un (s) = 0.895 >< 64^4(^) 1/2 ^ (59) 

is our final result for the positive signature mesonic singlet exchange in the t channel. 
The factor 0.895 indicates that the large A approximation has a 10% error in the 
normalization n(a s ). 



5 Conclusion 

In this paper we analyze the scattering amplitudes of quarks and anti-quarks through 
flavour exchange in the Regge limit. Both the normalization factors and the mesonic 
Reggeon trajectory are obtained. 

The negative and positive signatures of the fixed coupling constant case give nearly 
the same energy behavior at the SSC energies. The Regge trajectories are porpor- 
tional to Ja s (—t) which are non-linear in strong coupling constant. 

The inclusion of the effect of the running coupling constant reduces the magntiude 
of the amplitudes dramatically. The failure of small b expansion indicates that the 
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running coupling effect is not a small perturbation relative to the fixed coupling case. 
The change of the nature of singularty, from a square root branch cut to a simple pole, 
reflects the importance of the effect of the running coupling constant. Moreover, the 
position of the singularity shifts to the left and is a lot smaller than that of the fixed 
coupling case. Its smallness imposes a serious challange to observing experimentally 
the mesonic exchange in the parton level. 

Although in this paper we study just the positive signature in the running cou- 
pling case, as explained in Sec. 4, we do expect that positive and negative signature 
amplitudes have approximated normalization factors and trajectories. This enables 
us to make an estimation on the possibility of observing mesonic exchange in the 
parton level. 
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Figure Captions 

Fig. 1: Elastic quark anti-Quark annihilation. 

Fig. 2: The square-root branch cut of / + . C is the contour of the Mellin transfor- 
mation which lies to the right of the singularities. 

Fig. 3: The square-root branch cut of f . 

Fig. 4: Numerical solution (dot) of the Regge trajectory uor as function of strong 
coupling constant a s . The solid line is the numerical fitting given by eqs. (56, 
57). 

Fig. 5: Normalization factor n(a s ) from a s = 0.121 to 0.382 
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